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Abstract 

We study a monomer-dimer model with repulsive interactions be- 
tween the same species in one dimension. With infinitely strong in- 
teractions the model exhibits a continuous transition from a reactive 
phase to an inactive phase with two equivalent absorbing states. Static 
and dynamic Monte Carlo simulations show that the critical behavior 
at the transition is diff'erent from the conventional directed percola- 
tion universality class but is consistent with that of the models with 
the mass conservation of modulo 2. The values of static and dynamic 
critical exponents are compared with those of other models. We also 
show that the directed percolation universality class is recovered when 
a symmetry-breaking field is introduced. 



PACS numbers: 64.60.-i, 02.50.-r, 05.70.Ln, 82.65.Jv 

A monomer-dimer(MD) model was introduced by Ziff, Gulari, and Bar- 
shad to describe the oxidation of carbon monoxide on catalytic surface [0]. 
In this model, a monomer {CO) adsorbs onto a single vacant site, while a 
dimer (O2) adsorbs onto a pair of adjacent vacant sites and then immediately 
dissociates. A nearest neighbor of adsorbates, comprised of a dissociated O 
atom and an CO atom, reacts and forms a CO2 molecule and desorbs from 
the metal surface. In two dimensions, as the CO gas pressure is lowered, 
the system undergoes a first-order transition from a CO-saturated inactive 
phase into a reactive steady state and then a continuous transition into a 
02-saturated inactive phase. This continuous transition is shown to be in 
the same universality class as the directed percolation (DP) 0, |^, ^. In one 
dimension, there is only a first-order phase transition between two saturated 
phases. 



1 



Motivated by the monomer- dimer model, many related lattice models 
have been formed to study nonequilibrium phase transitions, e.g. the contact 
process(CP) the A model [Q, the interacting monomer-monomer model 
(IMM) d^, § , and so forth. The CP and A model are single-component models 
(there are two choices at a given site: vacant /occupied or healthy/diseased), 
while the MD and IMM model are multi-component models (three choices: 
vacant /monomer /dimer or the other monomer). A common feature of these 
models, whether single-component or mult i- component, is that they exhibit 
a phase transition from a reactive phase into an inactive phase of a single ab- 
sorbing state. The resulting critical behaviors are classified into the category 
of the so-called "DP conjecture" ^ which depicts that models exhibiting 
a continuous transition to a single absorbing state generally belong to the 
universality class of the directed percolation. 

The universality class for models with a single absorbing state is well es- 
tablished. But few studies have been made for models with more than one 
absorbing states. Recently Jensen and Dickman |jlO|, [ill] have extensively 
studied some nonequilibrium lattice models with infinitely many absorbing 
states, the pair contact process (PCP) and the reaction dimer (RD) model. 
Both models have a continuous transition from a reactive phase into an in- 
active phase with infinitely many absorbing states, which is shown, rather 
surprisingly, again in the DP universality class. 

This result might mislead the readers that the number of absorbing states 
is not relevant to the universality class of the absorbing phase transitions. 
This shows a sharp contrast to the case of equilibrium critical phenomena 
where the number of ground states plays a crucial role. The symmetry be- 
tween absorbing states may be more important than the number of absorb- 
ing states in determining the universality class. In the PCP, the infinitely 
many absorbing states are not equivalent probabilistically, i.e. some absorb- 
ing states can be reached more easily than other absorbing states by the 
PCP dynamics. So it is important to study a model with multiple equivalent 
absorbing states. 

A few models have appeared in the literature which have two equivalent 
absorbing states. Those are the model A and B of probabilistic cellular au- 
toma (PCA) introduced by Grassberger et al ||T^ , nonequilibrium kinetic 
Ising models with two different dynamics (NKI) |T5[, and the interact- 
ing monomer-dimer model (IMD) |T^. The PCA and NKI models are 
single-component models, while the IMD model is a multi-component model. 
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Critical behaviors of these models are different from DP but seem to belong 
to the same universality class (see Table 1). This result confirms our assertion 
that the symmetry between absorbing states is relevant to the universality 
class of absorbing phase transitions. 

Recently, the branching annihilating random walks (BAW) with offsprings 
have been studied intensively 0, |T9|, The BAW models exhibit 



a continuous absorbing phase transition from an active steady state into 
a single absorbing state. According to the foregoing DP conjecture, this 
phase transition should belong to the DP universality class. However, recent 
numerical investigations of the BAW models with even number of offsprings 
(BAWe) reveal that the BAWe models do not belong to the DP class but the 
same universality class as in the models with two equivalent absorbing states 
1|, m (Table 1). The BAW models with odd number of offsprings (BAWo) 



belong to the conventional DP class Dynamics of the BAWe models 

conserve the number of walkers modulo 2, while the BAWo models evolve 
without any conservation law. So the conservation law seems also relevant 
to the universality class of absorbing phase transitions. 

The remaining question is why the models with two equivalent absorbing 
states and the BAWe models with the mass (particle number) conservation of 
modulo 2 exhibit critical behaviors belonging to the same universality class. 
Grassberger described the PCA models in terms of kinks and showed that the 



number of kinks is conserved modulo 2 in the PCA models [|T^. The NKI 
model can be also described by kinks (or domain walls) and its dynamics 
conserve the number of kinks modulo 2 . The kink representation of the 
IMD model is complicated due to the presence of various kinks resulting from 
its multi-component nature. Each type of kinks has no conservation law but 
the total number of kinks is conserved modulo 2 again ||23| . So these models 
bear a close resemblance to each other, even though there is no exact mapping 
available between them. From the spirit of universality, one can argue that 
differences among the above models are just irrelevant details which do not 
affect the universal behavior, e.g. scaling exponents. 

In this paper, we take a careful look at these details and raise the question 
of whether models with the same conservation law in its kink representation 
always belong to the same universality class. This question is intimately 
related to the role of symmetry between absorbing states in determining the 
universality class of absorbing phase transitions. 

We consider the IMD model with infinitely strong repulsions between 
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the same species in one dimension. A monomer (A) cannot adsorb at a 
nearest-neighbor site of an aheady-occupied monomer (restricted vacancy) 
but adsorb at a free vacant site with no adjacent monomer-occupied sites at a 
rate k^. Similarly, a dimer {B2) cannot adsorb at a pair of restricted vacancies 
{B in nearest-neighbor sites) but adsorb at a pair of free vacancies at a 
rate ks- There are no nearest-neighbor restrictions in adsorbing particles of 
different species. Here we will consider only the adsorption-limited reactions. 
A nearest neighbor of the adsorbed A and B particles reacts immediately, 
forms the AB product, and desorbs the catalytic surface. Whenever there 
is an A adsorption attempt at a vacant site inbetweeen an adsorbed A and 
an adsorbed B, we allow the A to adsorb and react immediately with the 
neighboring B, thus forming an AB product and desorbing the surface. The 
system has no fully saturated phases of monomers or dimers, but instead 
two equivalent half-filled absorbing states. These states, I and II, comprise 
of only the monomers at the odd- or even-numbered lattice sites. A dimer 
needs a pair of adjacent vacancies to adsorb, so a state with alternating sites 
occupied by monomers can be identified with an absorbing state. 

In order to study the role of symmetry between these absorbing states, 
we introduce a symmetry breaking field h which favors one absorbing state 
over the other. This can be done by differentiating the adsorption rate of 
monomers at an odd-numbered vacant site and at an even-numbered one. If 
a monomer is chosen to adsorb on an even-numbered free vacant site, the 
adsorption-attempt is rejected with probability h {0 < h < 1). The case 
h = corresponds to the ordinary IMD model previously studied in details 
1^, |l^]. For finite h the monomers tend to adsorb more on an odd- numbered 



site than an even-numbered one. Therefore the absorbing state / can be 
reached more easily than the other absorbing state II by this dynamics. 

In this paper, we set h = 0.5 for convenience. Then the system can be 
characterized by one parameter p = k^/ (/ca+^b) of the monomer adsorption- 
attempt probability. The dimer adsorption-attempt probability is given by 
q = 1—p. We perform dynamic Monte Carlo simulations for this model. We 
start with a lattice occupied by monomers at all odd-numbered sites except 
at the central vacant site. Then the system evolves along the dynamic rules of 
the model. After one adsorption attempt on the average per lattice site (one 
Monte Carlo step), the time is incremented by one unit. 5000 independent 
runs are made up to 5000 time steps for various values of p near the critical 
probability pc- Most runs, however, stop earlier because the system gets 
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into the absorbing state I. We measure the survival probabihty P{t) (the 
probabihty that the system is still active at time t), the number of dimers 
N{t) averaged over all runs, and the mean-square distance of spreading B?{t) 
averaged over the surviving runs. At criticality, the values of these quantities 
scale algebraically in the long time limit 



p{t) ~ t-^ (1) 

N{t) ^ t\ (2) 
R\t) ~ t\ (3) 

and double-logarithmic plots of these values against time show straight lines. 
Off criticality, these plots show some curvatures. More precise estimates for 
the scaling exponents can be obtained by examining the local slopes of the 
curves. The effective exponent 5{t) is defined as 

_ logipwrnt/m 

log 

and similarly for rjit) and z{t). In Fig. 1, we plot the effective exponents 
against 1/t with h = 10. Off criticality these plots show upward or downward 
curvatures. From Fig. 1, we estimate pc — 0.414(1) which is much lower 



than in the ordinary IMD model without the symmetry breaking field [|T7 
(pc ~ 0.5325(5)) as expected. The scaling exponent is given by the intercept 
of the critical curve with the vertical axis. Our estimates for the dynamic 
scaling exponents are 

(5 = 0.160(5), 77 = 0.31(1), 2 = 1.40(20). (5) 

These values are completely different from the values in the IMD model 
without the symmetry breaking field (Table 1), but are in an excellent accord 
with those of the DP universality class; 6 = 0.1596(4), = 0.3137(10), and 
z = 1.2660(14) (see reference |^), although the exponent z shows very slow 
convergence. Similar results are obtained as we set the value of the symmetry 
breaking field h = 0.1. 

In summary, we have studied the IMD model with infinitely strong repul- 
sive nearest-neighbor interactions between the same species in one dimension. 
This system exhibits a continuous transition from a reactive phase into an 
inactive phase with two equivalent absorbing states. Its critical behavior is 
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different from the DP universality class but consistent with that of the PCA, 
NKI, and BAWe models. As soon as the symmetry breaking field is intro- 
duced, the two absorbing states of the IMD model become probabihstically 
not equivalent and this model exhibits a critical behavior in the DP univer- 
sality class. It implies that the symmetry between absorbing states plays 
a crucial role in determining the universality class. Our numerical results 
indicate that the system evolves mostly into one favored absorbing state, not 
bothered by the other unfavored absorbing state. So the critical behavior 
may be governed by a fixed point with a single absorbing state, which is 
associated with the DP universality class. This may explain why the PGP 
model with infinitely many absorbing states also belongs to the DP univer- 
sality class. For more numerical evidences, we are currently performing static 
Monte Carlo simulations for this symmetry-broken model and also studying 
the role of the symmetry breaking field in the PCA and NKI models. It is 
interesting to note that the total number of kinks in the IMD model is still 
conserved modulo 2 in the presence of the symmetry breaking field. So it is 
incorrect to say that all models with mass conservation of modulo 2 should 
belong to the same universahty class other than the DP class. 

This work is supported in part by NGN DIRECTED RESEARCH FUND, 
Korea Research Foundation and by the BSRI, Ministry of Education (Grant 
No. 95-2409). 
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Table Caption 



Table 1 : Numerical estimates of scaling exponents. r]d {r]i) is the exponent 

These estimates are taken from 



for the defect (interface) dynamics [ T7 
references H, |T|, |16|, |l 
errors in the last digits. 



\21^ . Numbers in parentheses represent the 



exponent 


PCA 


NKl 


IMD 


BAWe 


S 


0.283(16) 


0.27(2) 


0.29(2) 


0.285(2) 


Vd 






0.00(2) 


0.000(1) 


Vi 


0.272(12) 


0.30(2) 


0.285(20) 


0.282(4) 


z 


1.11(2) 


1.14(2) 


1.14(2) 


1.141(2) 




3.3(2) 


2.9(6) 


3.17(5) 


3.25(10) 




0.94(6) 


0.80(8) 


0.88(3) 


0.92(3) 
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Figure Captions 

Fig. 1 : Plots of the effective exponents against 1/t. Thrccc curves from top 
to bottom in each panel correspond io p — 0.413, 0.414, 0.415. Thick 
lines are critical lines (p = 0.414). 



9 



